PSR J0737— 3039B: A probe of radio pulsar emission heights 
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ABSTRACT 

In the double pulsar system PSR J0737— 3039A/B the strong wind produced by pulsar A dis- 
torts the magnetosphere of pulsar B. The influence of these distortions on the orbital-dependent 
emission properties of pulsar B can be used to determine the location of the coherent radio 
emission generation region in the pulsar magnetosphere. Using a model of the wind-distorted 
magnetosphere of pulsar B and the well defined geometrical parameters of the system, we deter- 
mine the minimum emission height to be '^ 20R]ys in the two bright orbital longitude regions. 
We can determine the maximum emission height by accounting for the amount of deflection of 
the polar field line with respect to the magnetic axis using the analytical magnetic reconnection 
model of Dungey and the semi-empirical numerical model of Tsyganenko. Both of these models 
estimate the maximum emission height to be ^ 2500 Rns- The minimum and maximum emission 
heights we calculate are consistent with those estimated for normal isolated pulsars. 

Subject headings: stars:neutron - pulsars 
1. Introduction 



The mechanism of pulsar radio emission and 
its origin within the pulsa r magnetosphere are not 
well understood fsee. e.g.. lHankins et al.ll2009( ). In 
general, it is thought to be due to coherent ra- 
diation from relativistic plasma streaming along 
open magnetic field lines. Radio emission height 
estimates can constrain the emission mechanism 
to some extent. In normal isolated pulsars, ra- 
dio emission heights have been estimated from 
their emission geometry inferred from radio polar- 
i zation combined with the rota ting vector model 
(|Radhakrishnan fc Cookell969l) and the pulse pro- 



come from well above the surface ( Rankinlll990t 
Mitra fc Rankinll2002l ). However, these techniques 



file w idths ("Gil fc Kiiak"l993'; 'Kiiak fc Gil'^ 
Kiiak fc Gil 2003). Gangadhara fc Gupta (2001| 



and lDvks et al.l (|2004l ) have also proposed a phase- 
shift method to determine the emission height. In 
general, these methods show that core component 
emission originates very close to the surface of 
the neutron star (NS), but the conal components 



are limited in that we observe only a small section 
of the magnetosphere of these isolated pulsars due 
to an unchanging line-of-sight. 

PSR J0737-3039A/B is a unique binary sys- 
tem which provides an excellent opportunity to 
study different emission regions due to relativis- 
tic spin precession, allowing us to observe differ- 
ent portions of the magnetosphere. We can also 
explore magnetospheric distortion, which affects 
the observed emission pattern. The two neutron 
stars of this system orbit each other in a 2.4-hr 
orbit; this is the only known pulsar binary system 
in which both ne utron stars have been detectable 
as ra dio pulsars (JBurgav et al.l 120031 : iLvne et al 



20041 ). The first-born recycled pulsar — hereafter 
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A — has a spin period of 23 ms and the second born 
pulsar — hereafter B — has a spin period of 2.8 s. 
The pulse profile of A has been stable since its 
discovery, but that of B has dramatically evolved 
through five years of observing, culminatin g in its 



radio disappearance in 2008 March (P erera et al. 
20101) . 



Due to the unstable features of B's pulse pro- 
files, both on long timescales and within a single 



orbit, it is challenging to understand the emis- 
sion geometry and the emission mechanism. By 
fitting a model to eclipses of A d ue to absorp- 
tion i n the magnetosphere of B, iBreton et al.l 
pOOSf ) constrained the geometrical parameters 
to be a = 70.9(4)° and = 130.0(4)° Q- Here 
a is the misalignment of the magnetic axis with 
respect to the spin axis and is the colatitude 
of the spin axis. Unfortunately, B shows very 
little radio polarization, making it impossible to 
constrain t he geometry from po la rization mea- 



surem ents ([Demorest et al.l 120041 ) . iBreton et al 



(|2008l ) also constrained the precessional phase of 
the spin axis, measured from our line-of-sight, to 
be (t)prec = 51.2(8)° at an epoch of 2006 May 2 
(i.e. MJD 53857) and found that this phase is 
changing with time at a rate of 4.8(7)° yr"-*^. This 
is consistent with the rate o f 5.061(2)° yr~^ pre- 



dicte d by gener al relativity (Barker fc O'Connelll 
19751 ). Recentlv. [Perera et al.l (|2010l) analvzed the 



pulse profile evolution of B and independently de- 
termined the above angles using a simple model 
ba sed on geodetic spin prec ession as proposed 
bv IClifton fc Weisberd (|2008[ ). According to this 
model, the beam must be elliptical and horse-shoe 
shaped in order to explain the observed single- to 
double-peak pulse profile evolution and the dis- 
appearance of radio emission. The estimations of 
the above angles i n this model are con sistent with 
those predicted bv lBreton et al.l ([20081), within the 
2-cr errors. From these studies, we believe that the 
geometrical parameters of B are well known. We 
can therefore use these values in this paper to 
determine the emission geometry. 

In this paper, we also explore the distortion of 
the magnetosphere of B. The almost edge-on or- 
bital pla ne of the system, w ith inclination angle 
of 88.7° ( Kramer et al.ll2006 ^. allows us to observe 
the eclipses of A with a duration of about 30 sec- 
onds. By considering the relative tr ansverse ve- 



lociti es of the two pulsars, 660 km/s (jLvne et al 



20041 ). the estimated size of B's magnetosphere 



is about 10% of its light cylinder radius of 
1.3 X 10^° cm . This implies that the wind of 
A compresses the magnetosphere of B and dis- 
turbs its polar cap (JLvutikovl 120041 ) . This is due 
to the small separation of the pulsars (^ 9 x 10^° 



^Here, and throughout the paper, the number in parentheses 
is the l-cr uncertainty in the last quoted digit. 



cm or 2.9 It-s) and the large spin-down lumi- 
nosity of A (5.8 X 10^'^ ergs~^) compared to B 
(1.6 X 10^° ergs~^). This is analogous to the dis- 
tortion of the Earth's magnetosphere due to the 
Solar wind. This interaction also results in an or- 
bital modulation of B; we detect bright emission 
from the pulsar only in two orbital phase regions of 
185° -235° (hereafter BPl) and 265° -305° (here- 
after BP2), and detect weak emission at phases of 
340° - 30° and 80° - 130° (Here, and through- 
out the paper, orbit a.1 phases a re me asured from 
the ascending node). iLvutikovl ()2005l ) claims that 
the pulsar has the same intrinsic radio intensity 
throughout the orbit and that the orbital modula- 
tion is due to the deflection of the magnetic polar 
field line with respect to the line-of-sight because 
of the influence of A. 

The wind interaction with the magnetosphere 
of B produces a bow shock between A and B; 
this is likely the boundary of the magnetosphere 
of B. The shape of this boundary depen ds on the 
orient ation of the magnetic axis of B. iLvutikov 
(|2004[ ) constrained the stand-off distance, or the 
distance from B to the vertex of the bow shock, 
to be 3.5 X 10^ cm if the bow shock interface is 
a perfect resistor and 4 x 10^ cm if it is partially 
resistive. These estimates inferred that the mag- 
netosphere of B is located deep within its light 
cylinder and the open and closed field lines have 
a more complicated structure than that of an iso- 
lated pulsar. Since the wind-interaction bound- 
ary model is very important to study the emission 
geometry of B, we derive it again in this paper 
with some improvements. This model describes 
the shape of the boundary for any orientation of 
the magnetic axis and allows us to model the open 
and closed field line structure more accurately. 

This results in a method to use the derived 
field line structure to estimate the radio emis- 
sion heights of B. Since the bow shock bound- 
ary is located deep inside the light cylinder, the 
correction due to rotation on the static dipole 
field, which introduced in retarded dipole field, is 
small (for retarded dipole field, refer appendix A of 
Dvks fc Harding 2004 ) . Therefore we assume non- 
rotating dipole field throughout the model. In our 
method, we assume that the emission comes from 
the direction tangential to the local field lines. We 
also assume that the emission comes from above 
the polar cap region, consistent with the narrow 



single- and double-peaked radio profiles. 

The distortion of the Earth's magnetosphere 
due to the solar wind has been studied using a 
large number of satellite observations and these 
data have been extensively modeled. Since the 
wind of A distorts B's magnetosphere in the same 
way that the solar wind does the Earth's, some 
models for the Earth's magnetic field line struc- 
ture can be used to study the distortion of the 
magnetosphere of B and to determine the r egions 
of radio emission. We use the lDungevI ( 1961 1 plan- 
etary m agnetosphere model and the iTsvganenkd 
( 2002ai a) Earth magnetosphere models to set an 
upper limit on the radio emission height of B. 

We present our observational data in section [5] 
and discuss observed mean pulse profiles. In sec- 
tion [3l we explain the boundary model which de- 
scribes the shape of the bow shock. Then we trace 
the dipole field lines and transform them from the 
co-rotating frame of the neutron star to the orbit- 
fixed frame. In order to derive the required angu- 
lar radius of the elliptical beam for the emission 
height calculation, we re-analyze the beam with a 
different geometrical framework in section 21 In 
section [31 we present the method which we use to 
estimate the minimum emission height and our re- 
sults. In section[6l we describe the maximum emis- 
sion height calculation using two different mag- 
netosphere models. Finally in section [71 we dis- 
cuss our results and compare them with predicted 
emission heights for other pulsars. We also discuss 
height estimations from other methods, conclud- 
ing that these are not applicable to pulsar B. 

2. Observations and pulse profiles 

We observed J0737 - 3039B with the 100 m 
Green Bank Telescope (GBT) in West Virginia 
since 2003 December 24 at multiple frequencies. 
Since 820 MHz is the most common, we use only 
those data in this analysis. This is the sam e 
data set which we reported in lPerera et al.l ( 20101 ). 
However, we use better time resolution pulse pro- 
files in this analysis compared to the previous pa- 
per. The data were taken using the GBT spec- 
trometer SPIGOT with sampling time of 81.92 /is 
until 2009 January. After 2009 January, the spec- 
trometer GUPPI was used with a sampling time of 
61.44 /j,s. All the data were dedispersed and folded 
using the pulsar analysis package SIGPROC, as- 



suming a dispers ion measure of 48.9 14 cm ^ pc 



( Lvne et al.l 120041 ) . The ephemeris of iLvne et al 



20041 ) was used until 2006 and sinc e then we have 



used the ephemeris of lKramer et al.l ( 20061) to form 
mean pulse profiles. 

The mean pulse profiles for BPl are shown in 
Figure m We aligned the peak of the profiles to the 
pulse phase of 0.5 at each epoch. Note that these 
pu lse profiles have a b etter resolution than those 



Perera et al 



(|2010t) . The second peak of the 
pulse profile can be hidden with the low time res- 
olution. Therefore, we use 1024 bins across the full 
pulse phase, resulting in an effective time resolu- 
ti on of 0.003 s in this analysis compared to 0.01 s 

(I20i3). 



Perera et al 



For example, the second 
peak of the pulse profile of MJD 53860 around 
pulse phase 0.5 2 in Figure [T] cannot be clearly seen 
in Figure 1 of IPerera et al.l ( 20101 ) on the same 
day with low time resolution. We use the pulse 
profiles from 23 days in section 21 including 16 
epochs in Figure [Tl in order to derive the beam 
shape. We include some low signal-to-noise data 
(e.g. MJD 53481) in this particular analysis be- 
cause the second peak became apparent around 
those days. Because these profiles appeared as 
single-peaked, we ignored them in the geometri- 
cal modeling of the previous lowe r time resolution 
study. As in IPerera et al.l ( 2010I ). we fit one and 
two Gaussians for each single and double-peaked 
pulse profile, respectively, and then calculate pro- 
file widths at different intensity levels. We use 
these data in section |4] to determine the geometry 
of the pulsar and the beam. 

3. Boundary model 

Due to the distortion of the magnetosphere, the 
properties of pulsar B are different from those of 
normal isolated pulsars. In isolated pulsars, we 
can determine the size of the magnetosphere by 
modeling the open and closed field lines, given the 
size of their light cylinder. However, as mentioned 
earlier, the magnetosphere of B is located deep 
inside the light cylinder and the structure of the 
open and closed field lines is more complicated due 
to the distortion from A's wind. 

In the first step, we approximate the struc- 
ture of the magnetosphere as a rotating vacuum 
dipole. Then we apply a simple mod el for the 
wind-magnetosphere interaction, as in iLvutikov 
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Fig. 1. — Mean pulse profiles for BPl on 16 different days including the very first observation which was 
made on 2003 December 24 (MJD 52997). All data have been observed at a frequency of 820 MHz. There 
arc 1024 bins across the entire pulse profile, resulting an effective time resolution of 0.003 s. Since predictions 
of absolute pulse phase are not available for these observations, we aligned the maximum peak to the pulse 
phase of 0.5 at each epoch. The horizontal and dotted lines show the baseline, or off-pulse mean, of the 
profile and the standard deviation of the off-peak region, respectively. The signal-to-noise ratio (S/N) of 
pulse profiles has decreased significantly from 284 (on MJD 52997) to 11 (on MJD 54552). 



(|2004l ). The wind of A creates a dynamic pressure 
on the magnetosphere of B. The magnetosphere of 
B creates a magnetic pressure which opposes the 
wind pressure of A. At some point, these two pres- 
sures equal each other; this interface is hkely the 
boundary of the magnetosphere of B. This bound- 
ary can be used to calculate the last open and 
closed field lines. We derive an expression for this 
boundary by equating the two pressures, 

#(f^)/(87r) = La cos" {-i)lA^ erf , (1) 



gion. This can be written as 



where f~^ is the distance vector of the boundary 
with respect to pulsar B, B{f^) is the magnetic 
field of pulsar B at I^b, La is the spin-down lu- 
minosity of A, rX is the distance vector of the 
boundary with respect to A, and 7 is the angle 
between the normal to the boundary and rX (see 
Figure [5]). The relative pressures lead to a bound- 
ary much closer to B. For that reason we simplify 
the problem by setting the distance of the bound- 
ary with respect to A equal to the distance be- 
tween the two pulsars and 7 equal to the angle 
between the normal to the boundary and the line 
connecting two pulsars. By assuming a magnetic 
dipole at the center of the coordinate system, we 
can write the magnetic field strength of the neu- 
tron star as B{r^) = {2>fB{fA ■ tb) — r^)/rB^, 
where r^ = m(cos S cos Qt, cos 5 sin fJi, sin 5) is the 
magnetic moment, n is the rotational angular fre- 
quency, and S is the angle between the magnetic 
moment and the line connecting the two pulsars 
(see Figure [2]). 

In isolated pulsars, we believe that the spin- 
down is caused by the power carried out along the 
magnetic field lines which are ope n with respect 
to th e light cylinder (Conto poulos fc: Spitkovsky 
20061 ). These open magnetic field lines start from 
the polar cap region of the neutron star surface 
and their shape changes with respe ct to the mag- 
netic inchnation. ISpitkovskvl (|2006l) proposed a re- 
alistic form of the spin-down luminosity of an iso- 
lated pulsar as a function of the magnetic inclina- 
tion angle. However, the open field line structure 
of pulsar B is somewhat different than for an iso- 
lated pulsar and, therefore, we need to define the 
magnetic fields lines with respect to the bow shock 
boundary. Therefore, we modified the spin-down 
luminosity equation that is given in SpitkovskyI 
( 20061 ) by including the area of the polar cap re- 



Lb = n^B^S^{l + sin^ a)/47r^c 



(2) 



where n is the rotational angular frequency, Bq is 
the polar magnetic field of pulsar B, 5 is the area 
of the polar cap region, and a is the magnetic in- 
clination. As we mentioned earlier, with previous 
geom etry models (JBreton et al.ll2008t IPerera et al 



2010t) . a is taken to be ~ 70°. Note_that_the origi- 
nal version of the equation (as in lSpitkovskvll2006 ) 
can be obtained by taking S to be the area of the 
polar cap of a dipole as defined by the open mag- 
netic field lines with respect to the light cylinder. 
In order to determine the boundary, we solve 
equation ([T]) numerically and then use equation 
([2]) to determine the value of the magnetic field. 
To do so, we simplified the problem to a 2D form 
in which the bow shock is represented by an equa- 
tion involving Xp and Zp and lies on that plane. 
The equation of the bow shock then has the form 
G{xp,Zp) = f{zp) — Xp, which must be solved in 
order to determine the shape of the bow shock. 
Let the radial vector tb = f{zp)xp + zip and the 
vector normal to the boundary n = V • G{xp, Zp). 
Then the dot product of these two gives the an- 
gle cos^(7) = 1/(1 -I- {df/dzp)"^) and this can be 
substituted in equation ([T|). According to our 2D 
form, we can write the magnetic moment m = 



I cos Sa 



I sin Szp and then derive the magnetic 



field of B at distance rs, B{fB), as a function of 
m, f{zp), z and S. Then equation ([T]) reduces to 
a first order differential equation of f{zp) and the 
solution determines the shape of the bow shock. 
First we assume an initial value for the magnetic 
moment, m, of B with a possible magnetic orien- 
tation, S, and solve the problem to determine the 
shape of the bow shock. Then we find which are 
the last closed field lines defined with respect to 
the bow shock and finally the shape and the area 
of the polar cap. 

We repeat the procedure for 16 values of the an- 
gle S between the magnetic moment and the line 
connecting the two pulsars, evenly spaced between 
and 7r/2, and we find the area of the polar cap for 
each of those orientations. In the estimation of the 
average area of the polar cap, we have weighted 
appropriately the fact that some values of S oc- 
cur more frequently than others during an orbital 
period. Using the spin-down luminosity given in 
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Fig. 2. — Geometry of ttie wind-magnetosphere interaction model. Pulsar A is located at a distance D 
(~ 9 X 10^° cm) away from B along the Xp-axis. The wind of A comes along the —Xp direction and is shocked 
near B. The physical interface has the shape of a bow shock and this is the boundary of the magnetosphere 
of B. fh is the position vector at any point on the interface with respect to B. r^ is the position vector of 
this point with respect to A and 7 is the angle between the normal to the boundary at this point and fX. 5 
is the angle between the magnetic moment axis rh and the line connecting the two pulsars. 



equation ([2]) with the timing-derived Lb, we find 
a new value for the magnetic field. We repeat this 
process with this new magnetic field until the value 
of the magnetic field converges. This happens af- 
ter five to ten iterations for an initial guess of the 
magnetic field within a couple of orders of magni- 
tude away from the convergence value. In order 
to represent the 3D version of the bow shock, we 
assume that it is axially symmetric around Xp. 

According to the best solution, the magnetic 
field of B is constrained to be Bb — 6.4 x 10^^ 
G, which is about a factor of two lower than the 
timing-derived value 1.2 x 10^^ G assuming a vac- 
uum dipole with a magnetic inclination of 90° 
(JLvne et al.l 120041 ). This new estimate is more re- 
alistic as it accounts for the boundary of the mag- 
netosphere as the bow shock and a realistic mag- 
netic inclination. The stand-off distance is con- 
strained to be either 3.8 x 10^ cm or 4.5 x 10^ cm 
for the cases when the magnetic axis is either nor- 
mal or parallel to the line connecting the two pul- 
sars, respectively. Thus, the size of the boundary 
depends on the orientation of the magnetic axis of 
pulsar B. Moreover, the shape of the bow shock 



depends on the orientation of the magnetic axis. 
The stand-off distance corresponds approximately 
to 1/3 of the light cylinder, thus for these dis- 
tances the relativistic modifications are minimal 
and do not change the value of the stand-off dis- 
tance by more than a few percent. For that rea- 
son we have chosen to calculate it using a vacuum 
dipole model rather than a more complicated ge- 
omet ry that takes in to account relativistic effects 
as in iDeutschI (|l955l ) . For simplicity, we assume 
the boundary is axially symmetric around the vec- 
tor connecting the two pulsars. A maximum devi- 
ation of roughly 20% of the actual shape from the 
symmetric case occurs when the angle 5 is 90°. 
Therefore, the shape of the boundary is sensitive 
to an angle of 5, having a range of [0°,90°]. If 
the angle 5 is greater than 90° , then the boundary 
considered the effective 5 of 180° — 5. For example, 
if i5 is 100°, then the effective 5 for the boundary 
shape is 80°. Therefore, with the assumption that 
the magnetic axis is nearly aligned with the line-of- 
sight at the radio emission detection, the effective 
angle 5 is small (^ 5° — 35°) in the orbital lon- 
gitude region of BP2. Thus, the deviation of the 



boundary model from the actual geometry of the 
boundary is small and the assumption of a sym- 
metric geometry is reasonable. However, the devi- 
ation in the orbital phase region of BPl is signifi- 
cant due to large effective angles of (5 (~ 35° — 85°). 
To derive an expression for the physical shape of 
the boundary, we examine different shapes which 
fit our results. As a preliminary fit, a parabola 
is a good guess, but a fourth-order polynomial 
describes the boundary better, yielding the mini- 
mum chi-squared value when we fit to our results. 
The coefficients of this polynomial describe the 
variation and are functions of the angle 6. The 
best-fit polynomial is 



x^ = aiS) + b{5){vl + 4) + c('^)(2^P + ^p)'' (3) 

with a B-centered coordinate system in which the 
Xp axis is towards A, the Zp axis is normal to the 
orbital plane, and the yp axis completes the right- 
handed coordinate system (see Figure [2]). The 
three axes have units in centimeters. The coef- 
ficients a(i5), b{5), and c{5) are 



a{5) = 

m = 

c{5) = 



0.83 - 0.01(5 - 0.06(52 _ 0.05(53 _^ o.03(5^ 



1.83 X 10 



-10 



-0.46 + 0.04(5 - 1.43(5^ + 1.96(5^ - 0.64(5^ 



5.45 X 109 

-0.48 - 0.03(5 + 2.15^2 _ 2.475^ ^ 0.74(5'* 
1.62 X 1029 



(f 



where the angle 5 is in radians and having a range 
of [0,7r/2]. Since they are functions of 5, the 
boundary changes slightly with spin and orbital 
motions, as well as over time due to precession. 

This boundary model is valid only up to 5 x 
10^ cm, or ~ 40% of the light cylinder radius, 
from B. Beyond this limit, the physical assump- 
tion of the dynamical pressure is incorrect because 
the wind pressure on the magnetic field should be 
zero when it is parallel to the boundary at large 
distances. Also we have assumed an undistorted 
magnetic field of B in the model and the distor- 
tions at large distances will be significant. 

In summary, we determined the boundary of 
pulsar B by assuming equilibrium between the dy- 
namical pressure of the wind of A and the mag- 
netic pressure of the field of B. The physical shape 



of the boundary is a bow shock and mathemati- 
cally we can represent it as a fourth-order polyno- 
mial. Moreover, this shape depends on the orien- 
tation of the magnetic moment axis with respect 
to pulsar A. Thus, the coefficients of the best fit 
polynomial depend on this orientation. 

3.1. Tracing the dipole field lines 

Due to the wind interaction with the magneto- 
sphere of B, it is complicated to understand the 
structure of the open and closed field lines. To 
determine the polar cap region that is required for 
the emission height estimation, we calculate the 
last closed field lines by tracing them referring to 
the derived boundary model. 

As is standard, we treat the magnetosphere of 
B as a magnetic dipole. For an isolated pulsar, the 
last closed magnetic field lines are defined as those 
that just touch the light cylinder and the ones in- 
terior to the last closed field lines are considered 
open field lines. In our case the boundary is not 
the light cylinder but the bow shock, with the last 
closed field lines defined as those that just touch 
this bow shock. The polar cap region is defined 
by these particular field lines and the shape of it 
can be determined by the locations where these 
field lines cross the neutron star surface. Defining 
the polar cap is important since we think that the 
radio emission is produced above this region. 

Unlike those of isolated pulsars, the magneto- 
sphere of B is not symmetric around the magnetic 
axis due to the shape of the boundary. This can be 
clearly seen by tracing the last closed field lines. 
In order to trace the field lines, we use the dipole 
field line equation in polar coordinates 



r = To sin (A) 

(j) = (j)Q 



(5) 



where r is the radial distance to a given point along 
the field line and tq is the field-line constant, or 
equatorial distance of the field line from the mag- 
netic axis. The angle A is the colatitude of a given 
point along the field line and 4>q is the azimuth an- 
gle, or the longitude of the given field line. Then 
the Cartesian components of a particular field line 
are written as 



X 

y 



r sin(A) cos{(j)) 
r sin(A) sin((/)) 
r cos(A) 



(6) 



where the z-axis of the coordinate system is 
ahgned with the magnetic moment axis and the 
other two axes are co-rotating with the neutron 
star. 

To include the misahgnment of the magnetic 
axis and also account for the spin phase, we trans- 
form equation ^ to another frame where the z- 
axis is aligned with the spin axis. In this frame 
(see Figure I31J a)), the Cartesian components are 



X cos a cos ( 



z sm a cos < 



X cos a sm < 



y smc 



+ y cos ( 



+ 



+ 



z sm a sm < 



spin 



(7) 



z cos a — X sm a 



where a is the angle between the magnetic axis 
and the spin axis and (j)spin is the spin phase. We 
measure the spin phase from the Xs axis, which 
means that it is zero when the magnetic axis is in 
the Xs-Zs plane. 

The spin axis is also associated with the co- 
latitude angle and the spin precessional phase. 
Geodetic spin precession changes the orientation 
of the spin axis with respect to our line-of-sight. 
In order to include the colatitude of the spin axis 
and effects of spin precession, we transform a par- 
ticular field line to another frame which is fixed 
with our line-of-sight. We choose the frame with 
z-axis parallel to the orbital angular momentum 
axis and a:-axis in the plane of the z-axis and the 
line-of-sight (see Figure[3jb)). The Cartesian com- 
ponents in this frame are 

Xo = XsCOsB cos (t>prec{t) - Vs siu 0prec (*) + 

Zs sin 9 cos 4>prec{t) 
Ho = a;^ COS 6* sin (jjprec (t) +ysCOS 4>prec (t) + 

Zs sin 9 sin (jjprecit) (8) 

Zo — Zs COS0 — XsSin9, 

where 9 is the angle between the spin axis and 
the orbital angular momentum axis and (j^precit) 



is the spin precession phase measured from the x- 
axis (i.e., 4>precit) = when the spin axis is in 
the plane of Xo-yo)- The spin precession phase is 
given by 



YpTecy") — ^'precl,^ -^0/ 



(9) 



where ^prec is the spin precession rate of B, which 
is 5.061(2)° yr~^, and Tq is the time when the 
</>prec(0 is zero, defined as the time when the spin 
axis is in the Xo-yo plane. By using the above set 
of equations, we can transform dipole field lines 
from the co-rotating frame of the neutron star to 
the orbit-fixed frame (a;o, j/o, -^o) where the Xq-j/o 
plane is in the orbital plane and Xo sin(z) is point- 
ing towards the observer, where i is the orbital 
inclination. 

In order to place the polynomial boundary in 
the {xo,yo,Zo) frame, we need to account for the 
orbital motion of B. Due to this motion, the ori- 
entation of the boundary changes with respect 
to our line of sight, because the location of A 
changes with respect to B. This relative motion 
changes the shape of the magnetosphere of B 
with respect to the line-of-sight and then the 
shape of the polar cap region. This results a 
variation in the emission height estimate across 
the orbit (more details are given in section [5]). 
We place the boundary model (Xp,yp,Zp) in this 
orbit- fixed coordinate frame with Zq \\ Zp and 
then rotate it corresponding to the orbital phase 

(florb with {Xp COS (florb - Vp sin (f)orb, Xp sin (j)orb + 

Up sin (j)orb,Zp), where (j)orb = V" + 90° (see Fig- 
ure ISlJb)). Here (j)orb is measured from the ascend- 
ing node and tp is the rotation angle between Xo 
and Xp axes. Note that this 90° angle is included 
to convert the rotation angle ip to orbital phase 
(porb as measured from the ascending node. Then 
we trace the last closed field lines according to the 
orientation of the boundary for the given orbital 
phase. For example, Figure 0] shows the confined 
magnetosphere in the boundary model with the 
last closed field lines on MJD 54050 (2006 Novem- 
ber 11). Here, we use our best-fit geometry of the 
pulsar from the beaming model that is described 
in sectionHl At this particular epoch, the spin pre- 



cession phase is • 



'prec 



46° and we use the orbital 



phase (j>orb = 223° and the spin phase (pspin = in 
the figure. Note that, at this spin and the orbital 
phases, the spin axis of the pulsar is in the Xp-Zp 
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Fig. 3. — Cartesian coordinate systems that transform magnetic field lines from the co-rotating frame of the 
neutron star to the orbit-fixed frame, (a) The spin phase 4'spin and the misalignment of the magnetic axis 



a are defined in the frame (xs,ys, Zg). The spin phase is measured from the x^-axis. Therefore 







is defined when the magnetic axis is in Xs-Zg plane, (b) The colatitude of the spin 9 and the precessional 
phase (fiprec is defined in the frame (xo, yo, Zo)- In this frame, the Xo-axis is in the plane of the Zo-axis and 
our line-of-sight (LOS). The spin precession (f>prec is measured from the Xo-axis. The frame of the bow shock 
boundary {xp,yp, Zp) is then placed in this orbit-fixed frame with Zp \\ Zq and rotate with an angle of ip 



defined from the Xo-axis. Then the orbital phase is defined 
node, io = 90° — i, where i is the orbital inclination. 



^orb 



— tp + 90° as it measures from the ascending 
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Fig. 4. — Confined rotating vacuum dipole in the boundary model on MJD 54050 (2006 November 11) — a 
view in the Xp-Zp plane (left) and a view in the Xp-yp plane (right). Pulsar B is located at the center of the 
coordinate system and the wind of A comes towards the —Xp direction. The dashed line shows the derived 
bow shock from the wind-magnetosphere interaction model and this models the open and closed field lines 
of the magnetosphere. The solid lines are the last closed field lines with respect to this boundary. The field 
lines that have higher latitude than these shown closed field lines are considered open field lines. The scale 
is in units of 10^ cm. Here, a = 61° and 9 = 138.5°; these are our best-fit geometrical parameters from 
section m This view corresponds to the orbital phase of 223°. Note, for clarity of plots, we take the spin 
phase as zero, so that the north pole of the pulsar is pointing below the Zp — plane. However, the shape 
of the magnetosphere changes with the spin and orbital motion and over time due to precession. 



plane. 

The spin of the magnetic axis also results in a 
change in the shape of the magnetosphere due to 
the misalignment of the magnetic moment. How- 
ever, the most important orientation of the mag- 
netic axis is when it reaches the closest approach 
to our line-of-sight (i.e. where we detect the emis- 
sion). In order to measure the point of this closest 
approach, or the impact parameter /3(i), we use 
the equations 



cos C(t) — sin 6 cos (j)prec {t) sin i + cos 9 cos i 



m = at) - a. 



(10) 



Here, ({t) is the angle between the spin axis and 
our line-of-sight at a given time and the other an- 
gles have the usual meaning. We calculate the spin 
phase which gives this particular closest approach, 
so that we can estimate the emission height only 
at this particular spin phase. 

Now we can transform field lines from the co- 
rotating frame of the neutron star to the orbit- 
fixed frame. By using the boundary model, we can 



trace the last closed field lines, which determine 
the shape of the magnetosphere. Due to different 
orientations of the magnetic axis with spin, the 
shape of the magnetosphere with respect to our 
line-of-sight changes, but we are interested only in 
the spin phase which gives the closest approach to 
us. Nevertheless, orbital motion and spin preces- 
sion change the shape and we need to account for 
these in emission height calculation. 

4. Re-analysis of the beam shape 

For the emission height calculation, we re- 
quire the an gular radius of the emission beam. 
Perera et al.l (2010 ') claimed th e beam shape of B 
to be e lliptical by modifying the lClifton fc Weisbere 
( 2008 ) geometrical framework. The model in 
Perera et al.l (|201C1I) used the 2D geometry of the 
beam after projecting it to a plane which is per- 
pendicular to the spin axis, resulting in a projected 
angular radius. We improve this model by using a 
more realistic 3D model in this paper in order to 
determine the actual angular radius. 

The angular radius of a circular emission beam 
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has been calculated using pulse p rofile widths an d 
an as sumed emission geo metry (JGil et alJ 119841 ). 
Since IPerera et al.l (J2010I ) showed that the shape 
of B's beam is not circular, we attempt to derive an 
approximate equation for an elliptica l beam in 3D 



and then follow the same analysis of iPerera et al 



pOlOh in order to determine the geometry. First, 
we construct the beam with a set of coaxial hol- 
low cones representing different intensity levels in 
a way such that the outermost one represents the 
lowest intensity level and then the intensity in- 
creases gradually inwards toward the maximum 
and then decreases until reaching the center of the 
beam. All these cones have cross sections with a 
constant ellipticity of aj/bj, where aj and bj are 
semi-major and minor axes, respectively, of each 
hollow-cone beam. Then, for any given longitudi- 
nal angular radius pij{t) (see Figure [5]), magnetic 
misalignment angle a, and impact parameter /3(i), 
the pulse profile width Wj{t) can be derived from 
spherical trigonometry (i.e. from spherical trian- 
gle FBD in Figure [5]) as follows 



, (t) — 2 arccos 



f cosjpi^jit)) - cos^ja + I3{t)) \ 



(11) 
where subscript j specifies different intensity lev- 
els of the pulse profile. Note that pij{t) is time 
dependent because the region where our line-of- 
sight cuts the beam is changing with time due to 
precession. In order to relate pi,j{t) with the ellip- 
tical beam shape, we derive an equation with the 
assumption that the cross-section of the beam is 
small enough to use ID trigonometry. Then the 
longitudinal angular radius pij{t) can be given as 
a function of (3{t), 



Pi,jW = -\/sin^(/Oaj) -cos2(paj)tan2(^(i)), 

(12) 
where paj is the angular radius across the semi- 
major axis of the beam (see Figure [5]) for a given 
intensity level and x — (^j/bj, which is a con- 
stant for all different intensity cones. This ex- 
pression shows that the minimum pij(t) of zero 
occurs when the line-of-sight just encounters the 
beam (i.e. /3(t) = Paj), resulting in Wj{t) = 0. 
The maximum pij{t) occurs when the line-of-sight 
crosses the center of the beam (i.e. f3{t) — 0), 



which leads to the maximum Wj{t). Therefore, by 
combining equation (fTTj) and p^ for a given a, 
/3(i), and Pa,j, we can calculate the pulse profile 
width Wj {t) for any given intensity level. 

In order to determine the geometry of B, we 
fit the model-predicted pulse profile widths to ob- 
served pulse profile widths of BPl at different in- 
tensity levels using th e same likelihood a nalysis 

(1201 d) . 



that was described in IPerera et al.l (|2010l ). The 
fit was done by searching the entire parameter 
space of a, 9, x, and Tq. For each combination 
of these parameters, we vary paj from 0° to 30° 
freely until we reach the best solution. Then we 
use a maximum likelihood analysis to determine 
the best-fit geometrical parameters. The best-fit 
model for BPl is shown in Figure |6l The es- 
timated geometrical parameters a — 6l.0°t-2Ao 
138.5°_/ ,10 are consistent with those 



and 9 
in 



.j^ 



Perera et all (|2010[ ) and iBreton et all (|2008h 



within the 2-a errors. The ratio x is constrained 
to be 2.6lQg, lower than the estimate of the pre- 
vious paper. Our new estimate is more believable 
because it has been derived from a full 3D view- 
ing model. In addition to these parameters, we 
derive Tq to be MJD 573991^5 (2016 January 12), 
which results in a precessional phase of 46° at an 
epoch of MJD 54050 (2006 November 11). This 
estimate is consiste nt with the value predicted by 



Breton et al.l (|2008l ) at the same epoch. However, 



this is a somewhat arbitrary parameter that can 
be chosen fro m our best-fit mode l. Note that the 
best-fit To in IPerera et ahl (|2010l) is about MJD 
33360 (1950 March 20), which results in a preces- 
sional phase of 73° at an epoch of MJD 54050. 
These two best-fit (pprec result in a shift of the 
hour-glass 2D pulse profile shape along the pre- 
cessional phase or time a xis (see Figure [5] o f this 
analysis and Figure 17 of lPerera et al.ll2010l) . 

The angular radius across the semi-major axis 
of the beam at the maximum intensity level Pa,ioo 
and 10% of the maximum Pa,io are constrained to 
be 9.9° and 14.3°, respectively. In order to de- 
termine the effective angular radius of the beam 
(more details are given in section [5|), we use these 
angular radii with the derived beam geometry. 

For emission height estimates for normal non- 
precessing pulsars, a beam shape is not essential 
because our line-of-sight always observes the same 
section of the emission beam. However, for pre- 
cessing pulsars, we must consider a beam shape 
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LOS trajectory 




LOS 



Fig. 5. — Elliptical emission cone in the frame (xs, j/s, Zg). The magnetic axis of the beam represents with 
rh and it is misaligned with the spin axis of an angle a. The trajectory of the line-of-sight across the beam 
due to rotation is denoted with DBE. The angle BOD is the longitudinal angular radius pij{t) of the beam 
for a given intensity level at a given time. The angle AOC is the angular radius across the semi-major axis 
of the beam paj for a given intensity level, which is time independent and fixed for the beam. The angle 
AOD is the effective angular radius of the beam Pe,j{t) for a given intensity level. 



in order to determine the emission height due to 
observing different sections of the emission beam. 
Therefore, in this particular case, we use our de- 
rived elliptical beam shape with the best-fit ge- 
ometrical parameters of B to estim ate emission 
height s in section [5] As we see in iPerera et al. 
(|2010l ). the pulse profile evolution is somewhat 
similar in both bright phases. Therefore, we use 
the above best- fit beam parameters from BPl 
in our emission height estimates for both bright 
phases. 

(|2010l) 



Perera et al 



reported that the radio 
emission of B disappeared in 2008 March, because 
the line-of-sight precessed away from the partially 
radio-filled elliptical beam. According to the same 
argument that the beam is not entirely radio loud, 
we can explain the radio disappearance with our 
new elliptical beam geometry. With the partially 
filled beam configuration, we can predict the reap- 
pearance of the radio emission as our line-of-sight 
precess back to the radio loud regi on of the beam. 
With the model described in Pere ra et al.l ([2010), 
the reappearance is predicted to occur in around 
2035 with the same part of the beam or in around 



2014 if the beam has two symmetric radio-filled 
portions. However, our new model describes the 
reappearance is predicted to occur in around 2024 
with the same part of the beam. Therefore, Fig- 
ure [6] changes with the partially-filled horse-shoe 
beam to non-detectable emission from 2008 to 
2024. If our new beam model is correct, then the 
beam should not have two symmetric radio-filled 
parts, because our line-of-sight crossed around the 
center of the beam when the disappearance oc- 
curred in 2008. Therefore, if there are two sym- 
metric parts, then we would be able to detect radio 
emission at present day. The two different predic- 
tions for the reappearance from the two models 
mainly occur due to two different best-fit Tq val- 
ues. These two different Tq values give two differ- 
ent solutions for the geometry of B, however, the 
best-fit values for parameters a and 9 are the same 
for the two models with the 2-ct uncertainty. 

5. Emission height calculation 

In order to estimate the radio emission heights 
of pulsar B, we use the previously defined bound- 
ary model, the field line tracing technique, and the 
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Fig. 6. — 2D pulse profile of the orbital longitude region 185° — 235° (BPl), assuming the elliptical hollow- 

I rr qO I c qO _i_n A 

cone beam. The best-fit geometrical parameters are a = 61.0°_2!4°: ^ — 138.5°_4;4o and x = 2.6lQg (errors 
are la). The corresponding Tq is MJD 57399 (2006 November 11), which is the time where the spin axis of 
the pulsar is in the plane of our line-of-sight and the orbital angular momentum ax is. Note that these a and 
9 values are consistent with previous results (jPerera et aLlbOlOt iBreton et al.ll2008l) . The dots are the widths 
at equal intensity levels of the observed pulse profile. Each horizontal row of dots represents an observation 
at a given epoch. Equal-intensity contours are the elliptical beam model-predicted pulse profile widths at 
different intensities. The intensity increases from the inner dashed line outwards until the first solid line, 
which is the intensity of the peak, and then decreases outwards again. The intensity levels are, from the 
inner dashed line, 80%, 90%, 100%, 90%, 80%, 70%, 60%, 50%, 40%, 30%, 20% and 10%. The vertical axis 
is calibrated in years and can also be considered the spin precession phase. Note that, this figure shows the 
emission from the full elliptical beam. If the beam is partially-filled according to 2008 radio disappearance, 
then the model predicts no radio pulse profiles from 2008 to 2024. 
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modeled geometry of the beam. As we mentioned 
earlier, we assume that the radio emission is pro- 
duced above the polar cap region and originates 
tangential to the local magnetic field lines. 

First, we need to determine the boundary of 
the polar cap region, given by the last closed field 
lines. This can be done by tracing the field lines 
with the derived boundary model as described in 
section 13.11 We assume that the radio emission 
comes from above the entire polar cap region, so 
that the outer edge of the pulse profile (i.e. 10% of 
the maximum intensity) corresponds to the region 
between the open and closed field lines, approx- 
imately equal to the last closed field line. Then 
we determine the emission height that originates 
from these last closed field lines. 

In order to determine the tangent to a partic- 
ular last closed field line at a given moment, we 
rewrite the coordinate transformations in section 



13.11 in matrix form (see, e.g.. iGangadharal 120041) . 
so that it is easy to evaluate the equations relevant 
for our calculation. First, we write the dipole field 
equation in the co-rotating frame of the neutron 
star as 



~^ cor = ro(sin'^ A cos (/>, sin'^ A sin </>, sin^ A cos A). 



(13) 



We then transform it to the orbit-fixed frame 



I nrh — A. ■ -D • T r 



(14) 



where. 



A = 



cos a cos ( 
cos a sin f 



^spin 



— sm ( 
COS(' 



'spin 
^spin 



sm a cos ( 
sin a sin < 



'spin 



'spin 



— sma 



COS a 



(15) 



and 

cos 9 cos (, 

B = I cos 9 sin q 

— sin I 



'prec 
'prec 



■ smf 



COS( 



'prec 
^prec 







sin 6* COS ( 

sin 9 sin c; 

cos 6 



^prec 
^prec 



(16) 

For the detection of radio emission, our line-of- 
sight must be parallel to the tangential vector of 
the given field line at a particular point. By locat- 
ing this point on the field line, we can determine 



the height of the radio emission. To evaluate the 
tangent to the field line, we take 'r't = dl^orb/dX. 
Then the unit vector along the tangential direction 
(ft — ~T^t/\^t\) can be written as. 



rt 



A-B- 



5-l-3cos(2A) 



3 cos (j) sin A cos A 
3 sin (h sin A cos A 



2 cos^ A — ; 



'A 



(17) 

The direction of the magnetic moment axis in the 
orbit-fixed frame can be written as 



■m^ A- B ■ . 



(18) 



Then we evaluate the angle between the direction 
of the magnetic moment axis and the vector tan- 
gential to the field line (r) at any time through 
the expression 



cos(r) — ft ■ TTT. 



l + 3cos(2A) 
v/l0 + 6cos(2A)' 



(19) 



This is the same as equation (8) in iGangadhara 
(J2004( ). At the point of detection of radio emission, 
we take the angle r to be equal to the efii'ective 
angular radius of the previously derived emission 
beam at the given time. According to our assump- 
tion that the outer edge of the pulse profile (10% 
maximum) comes from the last closed field line, 
we take r « pe,io(i), where Pe,w{t) is the efi^ec- 
tive angular radius of the beam (angle AOD of 
Figure [5|) at the 10% of the maximum intensity 
level corresponding to a particular impact param- 
eter I3{t). We can derive an equation for pe,io{t) 
by using the spherical triangle FAD of Figure [S] 
as 



cos{pe,io{t)) = cos(a) cos(a -I- /3(t)) + 

sin(a) sin(a + /3(t)) cos(wio(t)/2). (20) 

Here, wia{t), the model-estimated pulse width, 
can be found through equation (fTTj) and ([T2|) with 
the best-fit parameters a, 9, x, and Pa,io- The im- 
pact parameter /3(i) for the given time can be de- 
termined through equation ([9]) and ([T0| with the 
best- fit Tq. This wio{t) is simply the 10% pulse 
width of the 2D pulse profile given in Figure |6] at 
the given time. By simplifying equation (|19p. we 
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find an expression for A, which is the colatitude of 
the emission point. This expression can be written 

cos(2A) ^ -[cos(pe4o(i))y8 + cos2(pgio(i)) - 

sin2(pe,ioW)]. (21) 
This is same as equation (9) in iGangadharal 



(|2004[ ). so that the colatitude angle of the emission 
point in our complicated geometry is simplified to 
that of an isolated pulsar. Then the emission 
height of this point can be calculated by using the 
first equation of ([S]). However, determining the 
field line constant, rg, in this equation is difficult 
due t o the bow shock boundary and its varia- 
tion. iKiiak fc Gill (|l997l ) assumed that this tq is 
the light cylinder radius since the isolated pulsars 
that they have studied have low magnetic incli- 
nations. To determine rg for our particular case, 
we trace the last closed field line which is tangent 
to our line-of-sight at the closest approach of the 
magnetic moment with respect to the bow shock 
boundary. Then we use the dipole field equation 
to estimate the emission height. 

The orientation of the bow shock changes across 
the orbit with respect to our line-of-sight at the 
closest approach, so that the emission height 
changes with orbital longitude, because the last 
closed field lines are defined with respect to the 
bow shock. Moreover, when the central part of 
the radio beam crosses our line-of-sight, we will 
detect a double-peaked profile since the two edges, 
leading and trailing, of the beam cross our line- 
of-sight. Thus due to the different orientation 
of the magnetic moment axis at these two edges 
with respect to us, our line-of-sight is tangent to 
two different last closed field lines which have two 
different vq values. Thus the height of the emis- 
sion produced by the leading and trailing edges of 
the beam are different. This is shown in Figure [T] 
For example, the emission heights produced by the 
leading components of the beam are constrained to 
be in a range of [24±8,31±10] and [20±6,21±7] in 
NS radn (10 km) on MJD 54050 (2006 November 
11) for BPl and BP2, respectively. The heights 
of the trailing edge of the beam in BPl and BP2 
are constrained to be in a range of [15 ± 5,19 ± 6] 
and [21 ± 7,38 ± 12] in NS radn, respectively. The 
errors of the height estimates are calculated from 
the 1-CT uncertainties of the best-fit geometrical 



parameters from the beaming model. Thus, the 
uncertainty of the height estimate is in a range of 
[6,10] and [5,19] in NS radii for the leading and 
the trailing edge of the beam, respectively, across 
the orbit on this particular day. 

Also, due to precession of the spin axis, the 
emission height varies with time because the an- 
gle /3(i) varies with time. Again, there are two 
different heights for the leading and trailing edges 
of the beam. These are shown in Figure [5] and [3 
In BP2, the difference between the two heights is 
not as significant as in BPl because of the orbital 
position of B in BP2. In this region, pulsar B, A 
and our line-of-sight are roughly aligned, result- 
ing a less deviation in emission heights for the two 
edges of the beam. 

If the emission is produced from the boundary 
between the open and closed field lines, we can 
consider these estimates to be the actual emission 
heights for B. If the emission is produced elsewhere 
within the open field line region, these are lower 
limits on emission heights. 

6. Upper limit of the emission height 

We may also set an upper limit to the radio 
emission height by modeling the distortions of B's 
magnetosphere induced by the wind of A. Since 
both the magnetosphere and the wind are strongly 
magnetized, the distortions depend on the relative 
strengths of the magnetic fields and thus on the 
distance from the neutron star. Depending on the 
location of the radio emission region and the line 
of sight (and hence on the orbital position) an ob- 
server will detect different radiation signatures of 
the distorted magnetosphere. Inversely, by study- 
ing the orbital modulation and using a model of 
the distorted magnetosphere, we can deduce the 
location of the emission region. 

Similar to how the Sun distorts the Earth's 
magnetosphere, pulsar A produces a strong enough 
wind to interact with the magnetic field of B and 
shape its magnetosphere. The nature of this in- 
teraction will vary depending on the properties 
of the wind. For a wind with a substantial par- 
ticle fiux, the formation of a bow shock, simi- 
lar to case of the Earth, is expected. In this 
MHD confinement model, the shape of the Earth's 
magnetosphere is mostly determined by the pres- 
sure balance between the supersonic solar wind 
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Fig. 7. — The model-estimated radio emission height across the orbit on MJD 54050 (2006 November 11). 
The solid line represents the height of the emission produced by the leading edge of the beam and the dotted 
line represents that for the trailing edge of the beam. In order to consistent with observations, we used our 
best-fit geometrical parameters of the beaming model, a — 61° and 6 — 138.5°. At this epoch, (pprec is 46°, 
the impact parameter /3 is 3° and the corresponding spin phase is 126°. Note that the difference between 
the corresponding heights for the leading and trailing components of the beam is more significant in some 
parts of the orbit . The orbital longitude regions for BPl and BP2 are denoted with dashed and dot-dashed 
lines, respectively. 
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Fig. 8. — The model-estimated radio emission heights of the orbital longitude region of BPl vs time across 
one precessional cycle. This predicted variation of emission height with time is due spin precession, making 
different line-of-sight cuts across the radio beam. The solid line represents the emission height from the 
leading edge of the beam and the dotted line represents that for the trailing edge. Here, the orbital phase 
is fixed at 200° (BPl), but the spin phase changes with time corresponds to the /3 value. The time axis 
represents a full precession cycle, 71 yr. Here, we have used the same geometrical parameters as in Figure [7] 
and assumed a full elliptical beam, not a partially filled horse-shoe beam. If the beam is partially-filled 
then no radio emission is expected from 2008 to 2024. This is why the model still predicts radio emission 
at present-day MJDs. The emission height is zero (^ 2030 — 2067) when the line-of-sight is out of the radio 
beam. 
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Fig. 9. — As described in Figure 13 but for BP2 (orbital longitude is 290°). Note that in this bright phase, 
the emission height difference between the leading and trailing components of the beam is not as significant 
as in BPl. 



and the Earth's nearly dipolar field. This curved 
shape of the magnetosphere is reproduced well b 



current numerical models (jTsyganenkd 12002 



11 py 

m 



Tsvganenko k Sitnovll2007f ). 

On the other hand, for a strongly magnetized 
wind, reconnection between the wind and the com- 
panion's magnetic field lines must be considered. 
This results in an open structure for the whole 
magnetos phere, si n iilar t o the one originally pro- 
posed by iDungevI ( 19611 ) for planetary magneto- 
spheres. 

In the case of the double pulsar, it is unclear 
whether a MHD confinement model or a reconnec- 
tion model is more applicable due to the unknown 
composition of A's wind. However, we are mostly 
interested in the overall geometric structure of B's 
magnetosphere. For this purpose, it is sufficient to 
discuss magnetospheric structure in the most basic 
terms, relying on the models of planetary mag- 
netospheres. We consider two extreme, though 
complimentary, models of the Earth's magneto- 
sphere; the h ighly res i stive, analytical reconnec- 
tion model of iDungevI (|l961l l. hereafter D61, and 



after TS02. 

We use the D61 model of the open planetary 
magnetosphere as a simple analytical representa- 
tion of the distorted magnetosphere. The D61 
model states that the interplanetary magnetic field 
(IMF) may become reconnected with the terres- 
trial field along the day-side magnetopause, where 
the magnetopause is the boundary between the 
Earth's magnetosphere and the solar wind. This 
results in a distortion of the higher altitude regions 
of the inner magnetosphere. 

Forbes fc Speiseii (|l97l[ ) neglected the dynam- 



the fully screened, 3D numerica l hydrq dynamic 
confinement model of lTsyganenkol ( 2002ai bl). here- 



ics of the reconnection processes and modeled the 
Earth's magnetosphere as a linear superposition 
of two magnetic fields: the Earth's closed field 
and the solar wind's uniform field. Following this 
approach, we can represent B's magnetosphere as 
a simple addition of the pulsar's dipole field and 
the wind's uniform field. Similar to the IMF in 
the D61 model, we treat the magnetic field in A's 
wind as homogeneous in the vicinity of B, with 
the direction of the magnetic flux density vector 
perpendicular to the line connecting the two pul- 
sars (assuming a toroidal field). However, depend- 
ing on whether the large-scale toroidal field is pro- 
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grade or retrograde with respect to the orbit, the 
geometric structure of the magnetosphere can be 
significantly different (see Figure ITOl) . 

Alternatively, we can use TS02 for more pre- 
cise, three-dimensional modeling. This model is a 
data-based best-fit representation for the Earth's 
screened magnetosphere based on a large number 
of satellite observations. The model provides the 
option of adding the contributions from external 
magnetospheric sources such as the ring current, 
magnetotail current system, magnetopause cur- 
rents, and the large-scale system of field-aligned 
currents to the Earth's dipole field. 

We used the GEOPACK code repository devel- 
oped by Tsyganenko, with modifications to match 
the properties of the double pulsar system. In- 
stead of analyzing every current component in the 
TS02 model separately, we manipulated the global 
input parameters of the code which define the geo- 
metric structure of the magnetosphere. The shape 
and scale of the magnetosphere is controlled by the 
solar wind ram pressure and the dipole tilt only. 
Variations in the value of the ram pressure change 
the magnetosphere self-similarly. In the numerical 
model, the ram pressure is represented by the pa- 
rameter PARMOD(l) and has units in nPa. PAR- 
M0D(2) represents the disturbance storm time 
(Dst) index and has units in nT. The Dst index 
is a measure of the size and strength of the ring 
current, which contributes to the overall field con- 
figuration in the inner magnetosphere. The TS02 
model is designed in such a way that the structure 
of the magnetosphere within a stand-off distance 
from the Earth has a very small dependence on 
the components of the IMF. Hence, for simplicity 
we set the transverse components of the external 
field (PARM0D(3)=By and PARM0D(4)=S;,) to 
zero. 

We performed a visual fitting (see Figure [TT|) 
of the boundary produced by the TS02 code to 
the boundary produced by our theoretical model 
(equation ^ and Q in section [S]). We set PAR- 
M0D(3) and PARM0D(4) equal to zero and 
changed PARMOD(l) and PARM0D(2) until 
the shapes of the boundaries matched. The best 
fit values that we obtained are PARM0D(1)=8 
nPa for the solar wind ram pressure, PAR- 
MOD(2)=100 nT for the Dst index, and, by 
default, the zero transverse components of the 
IMF (PARMOD(3)=0 nT, PARMOD(4)=0 nT). 



This set of parameters produces a magnetosphere 
boundary with a stand-off distance of about 10.4 
stellar radii. In order to make the spatial scal- 
ing consistent with the properties of the double 
pulsar, we rescaled the stellar radius parameter 
RO from 1 to 0.0026. This change simply ensures 
that the stand-off distance is about 4000 stellar 
radii, which is the value assumed throughout this 
section. 

The obtained values of the parameters (PARM0D(1- 
4) and RQ) are not supposed to be physically real- 
istic; rather, they produce a magnetosphere with a 
shape and size (defined by the stand-off distance) 
that match the properties of the double pulsar. 
Moreover, there could be other successful fits since 
they are derived from the visual inspection of the 
boundaries (see Figure [TT|) . Nevertheless, using 
this particular set suits our purpose of modeling 
an approximate structure of pulsar B's distorted 
magnetosphere without using large computational 
resources. 

We employed the same criteria to estimate an 
upper limit for the emission height for both mod- 
els. We assumed that the elliptical emission beam 
is located close to the polar field lines, which are 
nearly aligned with the magnetic axis at r ^ Rlc, 
where Rlc is the light cylinder radius. The 
anisotropic distortion of the magnetosphere by the 
wind changes the location of the polar field line rel- 
ative to the undistorted magnetic axis (see Figure 
[T^ (b) and (d)). As a first approximation, the de- 
flection angle can be expressed as adefi ^ B^/Bp, 
where B^^ and Bp are the magnetic filed of pulsar 
A's wind and the magnetic field of pulsar B, re- 
spectively. Close to the neutron star's surface, the 
influence of the wind's magnetic held on the overall 
field structure is negligible. Therefore, adefi ^ 
at the surface and increases outwards as the wind's 
magnetic fleld becomes comparable to the pulsar 
field near the boundary. The ratio between the 
two flelds, and hence the amplitude of the deflec- 
tion, depends on the distance from the star as well 
as on S, the angle between B's magnetic axis and 
the line connecting the two pulsars. There is a 
certain height above which the distortion is strong 
enough to deflect the polar fleld line by more than 
the angular radius of the beam {pa,io — 14.3°), 
which is determined in section |4l Furthermore, 
if the component of the distortion perpendicular 
to the trajectory of the center of the beam in the 
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Fig. 10. — Geometric structures of Dungey-type magnetospheres in 2D. Magnetospheric models constructed 
by adding a wind's uniform field, in the direction shown by the arrows on the dashed lines, to a dipole field 
with northward orientation. Magnetic fields of 6.4 x 10^^ G and 10 G are assumed for the surface magnetic 
field of pulsar B and the wind magnetic field, respectively. The direction of the wind is the opposite of 
the arrow showing the direction to the pulsar A. Both (a) and (b) have the same direction of the wind 
and orientation of the dipole. Same is true for (c) and (d). However, two different possible orientations of 
the magnetic field in the wind result in very different overall magnetospheric structures, (a) shows smaller 
deflection of the polar field lines compared to (b). Magnetosphere in (c) is mostly open, whereas one in (d) 
is mostly closed. In the latter, a radius of the enclosed magnetosphere is about 4 x 10^ cm. 
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Fig. 11. — TS02 magnetosphere fitted to the theoretical boundary model. Field lines are plotted in solid. 
Out of which the green color is for the polar field lines. Dashed line represents the bow shock boundary 
model. Tilt of the dipole is 45° and Vg ^ 4 x 10^ cm. We fix the shape and scale of the boundary and fit 
the parameters of the TS02 model. 
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vicinity of our line-of-sight is large enough, then 
the emission beam can be pushed away from the 
Hne-of-sight to the extent that they do not inter- 
sect with each other for any spin phase (see Figure 
[T^ (b)). This wih render the emission beam un- 
seen. On the other hand, the opposite can be true 
if the distortion occurs mostly along the local tra- 
jectory of the beam. In this case, the visibility of 
the beam stays unchanged and a small shift in the 
spin phase, at which the emission beam is seen, 
might be the only observable imprint of the dis- 
tortions (Figure [T2| (d)). However, the geometry 
of pulsar B suggests that the former must be real- 
ized (Figure [I2l (a) and (b)). Therefore, in order 
to be able to detect pulsed radio emission from B, 
the deflection angle adefi should not exceed 14.3°. 
This places an upper limit on the height of the 
emission region. However, in order to be able to 
use this reasoning, a pulsar must be detectable 
through its pulsed radio emission. We therefore 
restrict our analysis to only BPl and BP2, the 
distinct radio- loud regions of the orbit. 

In order to make use of this criteria, we ana- 
lyzed the moments of the closest approach at the 
orbital phases within BPl and BP2. At the mo- 
ment of the closest approach, the magnetic axis 
is nearly aligned with the line-of-sight, making 
the angle 6 about the same as the angle between 
the line-of-sight and the line connecting the two 
pulsars. In turn, the latter is related to the or- 
bital phase as (90° — 0orb) (mod 180). Therefore, 
(5* ~ (90° - (t)orb) when -90° < (j)orb < 90° and 
<5* ~ {(jiorb - 90°) when 90° < cjiorb < 270°, where 
(5* is the value of 5 at the moment of the closest 
approach. For instance, at the orbital phase 185°, 
S^ - 95° and at 305°, 5, - 145°. Therefore, S^ 
varies within [95°, 145°] and [145°, 180°] for BPl 
and BP2 respectively. 

We calculated the deflection angle adefi using 
two different methods. We used a simple ana- 
lytical estimation for the simplified D61 model, 
while employing more complex numerical calcula- 
tions for the modified TS02 model. In the D61 
model, the system is characterized by three main 
parameters: the magnetic moment of the pulsar, 
the magnetic field in A's wind B^, and the angle 
between the two. Neither of the models depend on 
the absolute values of the pulsar and wind's mag- 
netic fields. In both cases, only the ratio of these 
two fields matter. This boils down to the assump- 



tion that the wind from A is strongly magnetized 
(i.e., the shape of the boundary and stand-off dis- 
tance is defined solely by the magnetic pressure 
balance). The pulsar wind is highly magnetized 
near the light cylinder. The particle component 
only takes over much further, closer to the termi- 
nation shock. Since the wind from A reaches B 
after only 1000 light cylinder radii, the assump- 
tion about its high magnetization is valid. There- 
fore, we can describe both fields by one dimen- 
sionless parameter: stand-off distance normalized 
to the stellar radius. At the moment of the clos- 
est approach, the magnetic moment of B is almost 
aligned with the line of sight. In turn, due to the 
peculiarity of the double pulsar, the line of sight 
is nearly parallel to the orbital plane. Thus, at 
the moment of the closest approach, B's magnetic 
axis and the magnetic field of A's wind are nearly 
coplanar. Therefore, for approximate estimates, 
the full three-dimensional analysis of the system 
is not necessary and we only carry out the calcu- 
lations for the two-dimensional configuration. 

In 2D polar coordinates (r. A) , the equation for 
the magnetic field lines reads as 



dr 



rdX 

TDtot 



(22) 



where A is the colatitude and is equal to 90° at 
the equator and to zero along the magnetic axis, 
which is the same colatitude angle that defined in 
equation ([5|) . The angle between the local tangent 
to the field line and the vector r can be approxi- 
mated as rdX/dr. In order to find the deflection 
angle of the polar field line due to the distortions 
by the wind, we consider the change in rdX/dr. 
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(23) 

As a superposition of the pulsar's dipolar and wind 
magnetic fields we take a simple addition of the 
two. Therefore, it follows from the equation (P^ 
that 



adefi 



Bl 



B\ 



BP + Bi; 



B? 



(24) 



Here, BP,BI and B:^ 



B^ are r and A compo- 



nents of the pulsar and wind magnetic fields, re- 
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Fig. 12. — Schematic view of the deflection of the emission direction by the wind, (a) and (c) show the 
configuration of the main axes of an isolated pulsar, with different orientation of the spin axis, however, (b) 
and (d) show the same picture with an addition of the wind from the companion, for the same configuration 
of axes as (a) and (c) respectively. Grey-shaded ribbon-type surfaces are the traces of the cross-section of 
the emission beam. If the line-of-sight intersects this surface then the observer detects the radio emission. 
In case of (a) and (b), the orientation of the line-of-sight with respect to the spin axis and magnetic axis is 
such that it does not intersect with the grey-shaded surface after distortion, as shown on (b) . This is not the 
case for (c) and (d) where the line-of-sight intersects with the grey-shaded surfaces, even after distortion. 
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spectively. The angle between the local compo- 
nents of the fields can be either (5* + A + 90° or 
90° — ((5* + A), depending on the orientation of the 
toroidal field. Thus, in the frame of the dipole, 
the magnetic field components can be expressed 
as follows 



/UoTO cos A 



r 


27r r3 


\^^J 


P _ 
A - 


/ioTO sin A 


(26) 


w 
r 


^TB^am{X + 5,) 


(27) 


w 
A - 


= TS^cos(A + (5,). 


(28) 



Here, to is a magnetic moment of pulsar B while 
Bq is a strength of the wind magnetic field. The 
wind magnetic field structure is believed to be 
toroidal. However, for the sake of simplicity, we 
assumed a locally uniform wind field across the 
whole magnetosphere of B, which is feasible since 
the radius of the light cylinder is much smaller 
than the orbital radius. In the equations ^7} and 
(PS)). B^ and B™ can switch signs depending on 
whether the large scale toroidal field of the wind 
is prograde with respect to the orbital motion of 
the pulsars or retrograde. Mathematically it is 
equivalent to replacing 5, with 5, + 180°. Below, 
we derive the approximate expression for adefi for 
the upper signs in the equations ((27)) and (|28)) and 
only in the end substitute S^, + 180° instead of (5, 
to account for both cases. 

It is reasonable for our case [adefi < 14.3°) to 
limit our estimations to the field lines close to the 
magnetic axis, i.e. A ^ 0°. Then we can rewrite 
the equations (|25 |) - (f28| in the following way 



1 



27r r3 


(29) 


Bl^Q 


(30) 


B^!' = -B^ sinS, 


(31) 


B^ = -B^cosS,. 


(32) 



Here, we are only with left the upper signs in the 
second pair of equations. We get an approximate 
expression for the defiection angle by substituting 
(El-lHl into (El 



Odefl = 



B^cosS. 

Mo™ 



27r r^ 



B^sinS^ 



(33) 



We can rewrite equation (I33p in terms of normal- 
ized distance r = r/rg, where Vg — (/ioTO,/27ri?^) 
is a stand-off distance (r^ ~ 4 x 10^ cm as esti- 
mated from our boundary model). 



Oidefl 



r cos i5* 
1 -I- f^ sin (5* ' 



(34) 



From the criteria for the pulsed radio emission de- 
tectability, it follows that the values of r for which 
C(defi > 14.3° must be excluded as possible emis- 
sion heights. To find such values of f for any (5*, 
we use the condition that the absolute value of the 
right-hand side of the equation (j34p must exceed 
14.3°. Hence, in the case of a prograde toroidal 
field we have 



1 



r^ sin (5* 



> 14.3°. 



(35) 



By replacing (5* with S^, + 180°, we get the de- 
tectability criteria for the retrograde configuration 



-r^ cos (5* 



1 — r^ sin (5» 



> 14.3°. 



(36) 



The minimum of the values of f that satisfy equa- 
tions (|35|) or (|36|) for any S^ corresponding to BPl 
and BP2 is the best upper limit we can put on 
the emission height with this approach. In Fig- 
ure [131 we show the solutions of the equation ([35]) 
and (p6l) , represented by the shaded areas over the 
contours of constant deflection angle of 14.3°. As 
we can see in Figure [T31 both orientations of the 
wind magnetic field produce almost the same up- 
per limit, which is ^ 2500-Rjvs- This value cor- 
responds to (5* ^ 170° — 180°, i.e. when pul- 
sar B, pulsar A and the Earth are nearly aligned, 
where pulsar B is in BP2 region. This is consis- 
tent with the method by which we estimate an 
upper limit. At superior conjunction, the wind 
magnetic field is perpendicular to the magnetic 
axis, resulting in the largest distortion of the polar 
field lines. Therefore, the distance from the star 
surface, above which the deflection is more than 
14.3°, is smallest at superior conjunction. 

We can use the same criteria to set an upper 
limit on the emission height using the numeri- 
cal TS02 model. We require that the numerically 
calculated distortion angle of the polar field line 
must not exceed the angular radius of the beam 
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Fig. 13. — Excluded values of the radio emission height calculated by using the analytical D61 model. 
The curtain-shaped shaded areas represent the domains for emission height and S* for which the deflection 
angle exceeds 14.3°. The area with blue stripes leaned to the left corresponds to the prograde configuration 
(equation[35]), while the area with green stripes leaned to the right corresponds to the retrograde configuration 
(equation [36]) ■ The shaded areas between the vertical dashed lines represent the values of S* valid for the 
analysis (i.e. those corresponding to BPl and BP2). The black triangle marks the minimum value of 
r^ 2500NS radii within the shaded range of S*. Thus, 2500NS radii is the best upper limit for the emission 
height within the D61 model. 
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{pa. 10 — 14.3°). This allows us to find the maxi- 
mum emission height for each value of J, . 

Using the modified TS02 code, we can trace any 
field line (particularly polar field lines) of B's dis- 
torted magnetosphere for any orientation of the 
magnetic axis with respect to A's wind. This 
means that we can calculate the deflection angle 
for any (5* at any altitude. The 14.3 contour on 
Figure [T3] shows the altitudes at which the deflec- 
tion angle equals 14.3° for all values of (5*. The al- 
titude with the lowest value amongst others is the 
best upper limit we can put on the emission height. 
Since we can only consider ^*'s corresponding to 
BPl and BP2, the resultant best upper limit of the 
emission height would be 2500i?jvs for S^ ^ 95° 
(orbital phase of 185°, which is in BPl)(see Figure 
[II. 

In summary, we adapted the models of the 
Earth's distorted magnetosphere to the double 
pulsar system, based on the similarities between 
the two. We adjusted the spin-orbital and mag- 
netic field parameters corresponding to the obser- 
vational data. Both magnetospheric models, the 
analytic D61 and numerical TS02, draw simplified 
and very extreme pictures of the double pulsar. 
Nevertheless, both models offer an improvement 
over a simple dipole and allow us to set an upper 
limit on the altitude of the emission region using 
a novel technique. Moreover, they would greatly 
compliment each other if somehow unified into one 
model. 

We used the criteria of the pulsar emission de- 
tectability to estimate an upper limit of the emis- 
sion height. This requires the distortion angle of 
the polar field line not to exceed the angular ra- 
dius of the beam, derived from the observational 
data. 

Moreover, we arrived at similar results by em- 
ploying two very different models; the analyti- 
cal, highly resistive Dungey type model and the 
numerical, fully screened modified TS02 model. 
Both approaches led to the conclusion that B's 
radio emission is generated within the inner 22% 
of the light cylinder. 

7. Discussion 

The determination of radio pulsar emission 
heights is important for understanding their emis- 
sion mechanisms. Pulsar B of the double pulsar 



system provides a unique opportunity to study 
different emission regions of the magnetosphere 
due to precession. Also, the magnetosphere is dis- 
torted, exhibiting a complicated field line struc- 
ture, due to the wind of A. These distortions de- 
pend on the orbital and rotational phases of B. 
Observations of these distortions, not observed 
in isolated pulsars, via the orbital variations of 
the radio intensity of B allow us to pinpoint the 
location of radio emission. 

We have applied a simple wind-magnetosphere 
interaction model to determine the boundary of 
the magnetosphere of B. The best solution de- 
scribes the shape of the boundary as a polynomial, 
with coefficients dependent on the angle between 
the magnetic axis and the line connecting the two 
pulsars. Furthermore, the boundary was not ax- 
ially symmetric, but for simplicity we modified it 
to be symmetric. The stand-off distance ranges 
from 3.8 — 4.5 x 10® cm according to the orienta- 
tion of the magnetic axis with respect to the line 
connecting two pulsars, resulting in a size of the 
polynomial boundary of less than 30% of the light 
cylinder radius. However, this size is three times 
larger than the size inferred from eclipses of A. 
Thus, a possible explanation for this is that the 
particle density in the magnetosphere of B falls 
off significantly as a function of the radial distance 
from the center of the pulsar, so that the radiation 
of A penetrates the outer regions of B's magneto- 
sphere. 

Moreover, the variation of the boundary will 
change the shape of the open field line region. As 
a result, the spin-down luminosity of B can vary 
slightly due to the variation in the area of the polar 
cap. This causes a 1.5% periodic variation in the 
spin-down luminosity. It can also lead to a cor- 
rection on the spin phase , but this is very small 



(jGourgouliatos et al.ll2011[ ). Thus we did not con- 
sider this effect in our model. 

As we determined, the range of the allowed 
emission height depends on the orbital motion 
due to the relative orientation of the magnetic 
axis with respect to the boundary. Also, preces- 
sion changes the location of the spin axis, so that 
the emission height changes with time. In both 
of these variations, we have been calculating the 
emission heights for both the leading and the trail- 
ing edges of the conal elliptical beam. For a nor- 
mal pulsar with its light-cylinder boundary, these 
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Fig. 14. — Fixed value contours for the deflection angle calculated by using the modified TS02 model. Each 
contour shows the upper limits to the emission region height for the corresponding value of the deflection 
angle and S*. The shaded area between the vertical dashed lines represents the values of 5* valid for the 
analysis (i.e. those corresponding to BPl and BP2). The black triangle on the 14.3° contour marks the 
minimum value of 2500 NS radii within the shaded range of 5*. Thus, 2500 NS radii is the best upper limit 
for the emission height within the modified TS02 model. 
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two edges give the same height due to cyhndrical 
symmetry. When the impact param.eter is equal 
to the angular radius of the beam across the semi- 
major axis, we would detect a single-peak profile, 
resulting in one emission height. Figure [7] shows 
that the relative heights of emission due to the 
two components of the beam switch in the two 
bright phase regions since the orientation of the 
boundary changes throughout the orbit. These es- 
timated emission heights are about 1% of the light 
cylinder radius or 4% of the stand-off distance. 
Again, these should be considered lower limits if 
the emission does not originate on the last closed 
field line. Moreover, the analytical and numerical 
approaches to the upper limit estimate lead to the 
conclusion that pulsar B's radio emission is gener- 
ated within 22% of the light cylinder. 

In normal pulsars, radio emission heights have 
been calculated by us ing their geometry and the 
pulse profile widths (|Kiiak fc Gi]||l997t ). These 
range from about 10 to 100 Rns, less than 10% 
of the light cylinder radius. Our emission height 
estimates are consistent with these results. Thus 
the radio emission produced by B likely has the 
same mechanis m as for i s olated pulsars, which is 
consistent with iLvutikovl ( 20051 ). Most theories of 
pulsar radio emission place the generation region 
close to the star, typ ically within one stellar radius 
(e.g. iMelrosd 119951 ). In contrast, a model based 
on the anomalous cyclotron-Cherenkov resonance 



(JMachabeU fc Usovlll979l:lLvutikov et al.lll999l) re- 
quires emission to be generated much higher up in 
the magnetosphere, at hundreds of stellar radii. 
The fairly high emission altitudes of radio emis- 
sion inferred in the present paper are consistent 
with the latter models. 

The magnetospheres of pulsar s can be distorted 
due t o rotation as proposed in lOvks fc Hardina 
(J2004l ) , resulting in a rotational sweepback of the 
magnetic field lines. They found that at low al- 
titude the rotation deflects the local direction of 
the magnetic field line by at most an angle of the 
order of (r / R^c)^ j where r is the radial distance 
of the field line. We applied this rotational sweep- 
back model to pulsar B along with our boundary 
model and found that the rotational sweepback is 
very small, because the deflection of the magnetic 
field line from its local direction is of order 0.1 
radians. The estimated upper limit for the emis- 
sion height implies that this effect is less than 0.05 



radians, negligible compared to the distortions by 
the wind. However, this effect is significant when 
the radial distance of the field is close to the light 
cylinder, so that it is important in normal pulsars. 

The relativistic phase-shift method can be used 
to determine th e radio emission height s of pu lsars 
as describe d in [ Gangadhara fc Guptal (J2001I ) and 
Dvks et al.l (|2004[ ). This method uses aberration 



and retardation effects to explain the observed 
pulse phase shift of pulse profiles containing core 
and conal components. We applied this method to 
the double-peaked pulse profiles of pulsar B to es- 
timate the emission height. Because pulsar B has 
only a conal component, we assumed that pulse 
phase of zero was at the minimum between the 
two peaks. Then the phase-shift is measured from 
the two peaks, leading and trailing, with respect to 
this reference phase. The calculation shows that 
the phase-shift method does not work for pulsar B. 
For example, the emission height on MJDs 53860 
and 53939 is 6 and 23 NS radii, respectively. On 
MJD 54050 it is zero due to zero phase shift. Also 
on MJD 54400, the phase of trailing component is 
larger than the absolute phase of leading compo- 
nent, so that the emission height becomes nega- 
tive. The reason of these calculated height fluctu- 
ations is that the pulse profile of B is not stable 
and varies significantly. Therefore, it is difficult to 
measure the shift in pulse phase accurately. Also, 
as there are only two peaks in the pulse profile 
of B, the determination of the pulse phase zero 
reference point is difficult. Therefore, the mea- 
sured shifts and then the emission heights may 
not be correct, concluding that this method can- 
not be used to constrain emission heights of pulsar 
B. However, this is a useful method of estimating 
emission heights of normal pulsars which have sta- 
ble pulse profiles with both core and conal em is- 
sion components ( Gupta fc Gangadharall2003h . 

In summary, by using the method presented in 
this paper, we can place limits on the radio emis- 
sion height for any pulsar with well-determined 
emission geometry. The advantage of this method 
is that by estimating the field-line constant by 
tracing the magnetic field lines, we can constrain 
the emission heights of pulsars which have high 
magnetic inclinations that make them unsuitable 
for the other methods. Our radio emission height 
estimations for pulsar B will be useful for future 
studies and in particular can be used to constrain 
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propo sed g eometrical models such as iLvutikov 
(|2005l ) and iFreire et all (|2009l ) in order to accu- 
rately explain the observations. 
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